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Abstract 



We study the oscillatory flux dependence of the supercurrent in a thin 
superconducting loop, closed by a Josephson junction. Quantum fluctuations 
of the order parameter in the loop affect the shape and renormalize the am- 
plitude of the supercurrent oscillations. In a short loop, the amplitude of the 
sinusoidal flux dependence is suppressed. In a large loop, the supercurrent 
shows a saw-tooth dependence on flux in the classical limit. Quantum fluc- 
tuations not only suppress the amplitude of the oscillations, but also smear 
the cusps of the saw-tooth dependence. The oscillations approach a sinusoidal 
form with increasing fluctuation strength. At any finite length of the loop, the 
renormalized current amplitude is finite. This amplitude shows a power-law 
dependence on the junction conductance, with an exponent depending on the 
low-frequency impedance of the loop. 

PACS numbers: 74.50 +r, 74.40.+k, 74.20.Fg 
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I. INTRODUCTION 



Quantum effects in ultrasmall Josephson junctions have been studied intensely now for 
more than a decade, both experimentally QXJ] and theoretically [0]. The most important 
manifestation of quantum fluctuations is the well-known macroscopic quantum tunneling of 
the phase across a current-biased junction [p]]. This phenomenon leads to the observation of 
events of quantum phase slip at a bias that is relatively close to the critical current. 

More generally, macroscopic quantum tunneling causes a finite voltage to appear at 
any finite current. The corresponding I — V characteristic may be nonlinear, and in the 
limit of zero temperature one finds V = AJ 7 , with an exponent 7 which depends on the 
impedance of the leads. This relation may be derived in the framework of the Caldeira- 
Leggett model || that treats quantum transitions between neighboring minima of a tilted 
washboard potential in the presence of a dissipative "environment". It has been shown in 
Refs. that the coefficient A is proportional to the square of the "bare" (i.e., unaffected 
by the environment) tunnel matrix element for transitions between the two minima of the 
potential. The dual result for a voltage-biased junction shows that the DC current in 
such a junction is proportional to the square of the Josephson energy of the junction, also 
unrenormalized by the environment 

Actually, these two complementary results follow from very similar treatments of two 
closely related models, the Caldeira-Leggett model ||, and the electromagnetic environ- 
ment model P,pH|: The effective boundary conditions for the quantum fluctuations of the 
"environment modes" in these treatments do not depend on the Josephson energy of the 
junction itself. This approach (adequate for most existing experiments) is absolutely le- 
gitimate in the case of weak fluctuations of the phase of the order parameter across the 
junction. A more cautious analysis is needed, though, if these fluctuations are strong. In- 
deed, quantum phase fluctuations in a one-dimensional (ID) superconductor are known to 
diverge logarithmically with length. If these fluctuations would result in a diverging random 
phase across the junction, the Josephson energy of the system would average to zero, and 
neither of the quoted results would be true. 

We will show that a finite renormalized Josephson energy can arise because the junction 
itself affects the fluctuations of the environment [pjj . Simultaneously, the modes of the 
environment renormalize the plasmon oscillations in the junction. This mutual influence of 
different parts of the circuit makes the separation of it on two entities - the junction with a 
fixed capacitance, and the environment - to be somewhat a matter of convention. 

We consider a thin superconducting loop which contains a Josephson junction. The 
"environment modes" of the loop consist of propagating plasmon modes with a soundlike 
dispersion |I2|]. The most straightforward way to observe a possible renormalization of the 
Josephson energy, is to phase-bias the junction by threading a flux $ through the loop and 
measure the loop magnetization, which is proportional to the Josephson current J = J($). 
The renormalization leads to values of the critical current J c which are smaller than one 
would expect from the mean-field result, J° = 7rAG/(2e), where G is the conductance of 
the junction, and A is the superconducting gap in the loop. For a large loop, it also changes 
the flux dependence </($)• 

The paper is organized as follows. The model for the loop with the junction is presented 
in Section ||, followed by a qualitative discussion of the dependence of J on $ in Section |TJ. 
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The problem of the renormalization of the Josephson coupling is treated in Section [TV] where 
we make use of the similarity to the problem of quantum Brownian motion in a periodic 
potential El. The effect of the macroscopic quantum tunneling on the phase- dependence 
of the Josephson current for a relatively large loop is discussed in Section [V]. We employ the 
analogy with the problem of pinning of a ID charge density wave jrjj], or Wigner crystal [jl5 



and use an instanton approach |I| to calculate </($)• Concluding remarks can be found in 
Section [VI|. 



II. THE MODEL 

We consider a superconducting wire of length L and small cross-sectional area S = ax a, 
which is embedded in a medium with dielectric constant e. The wire is closed by a Josephson 
junction to form a loop. The bare Josephson energy of the junction is Ej = 7iAhG/(4e 2 ), 
and its charging energy is E c = 4e 2 /C. Perpendicular to the loop, a magnetic field 7i is 
applied, such that a flux T-LL 2 /{An) threads the loop. We introduce the corresponding phase 
$ = HL 2 / (2$ ), where <3>o is the superconducting flux quantum. 

The low energy excitation spectrum of this system can be described in terms of the phase 
of the order parameter <p(x) by the Lagrangian £ = K — U, where 

K 2 r,.,-ua , R 2 im~ ¥>«>)] 



Jo 2e c LVA n 2E C K ' 



and 



Jo 8m \ ox L J 

-E jCos[cp(L)-ip(0)}. (2) 

Here, l/e c = [8e 2 ln(R/a)/e}~ 1 is the characteristic inverse charging energy per unit length of 
the loop (R is the distance to a metallic screen fL7|), n s is the density of the superconducting 
condensate, and m is the electron mass. 

The first term on the right hand side of Eq. ([!]) together with the first term on the right 
hand side of Eq. (0) describe the propagating plasma mode ]12| along the loop ||18|| . They 



correspond to the electrostatic energy stored in the plasmons and to the energy associated 
with the supercurrent in the wire, respectively. In the latter energy, we included only the 
kinetic inductance, assuming the wire is thin and electrodynamic effects are weak. The 
plasma mode is characterized by a linear dispersion relation u>(k) = v p \k between frequency 
u(k) and wave vector k, where the plasma velocity is given by 



' e c n s S _ _a_ 2 \n(R/a) 
Am ~ C 2\ L \ ns U 



Here, c is the speed of light and \l = \Jmc 2 / (47m s e 2 ) is the London penetration depth of 
the wire. At temperatures T much smaller than the superconducting gap A, the plasma 
mode involves oscillations of the supercurrent only, and damping due to thermally excited 
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quasiparticles is negligible. Retardation effects were neglected in the derivation of (|l|) - @, 
i.e., we assume that v p \ <C c; hence we require a -C A^. 

The remaining terms in (|1|) and (^|) refer to the Josephson junction. For simplicity, 
we will completely neglect the junction capacitance, C = 0, throughout this paper. This 
approximation corresponds to neglecting the last term in comparison with the first one in ([!]) 
for all the relevant scales of the phase variation. The shortest scale in the time dependence of 
ip is h/A, and correspondingly the smallest part of the ring involved in the phase fluctuation 
is hv p \/A. Therefore, we may set C to zero if 

Ec mv pl 

A {Hn s Sy y ' 

An important length scale in our model is determined by the length L* at which the 
energy of supercurrents in the loop and the Josephson energy of the junction are of the 
same order. It is given by 



_ h 2 n s S _ hgv pl 



where g is defined through 



1 = 4mvpi_ = q ^Xl / 21n (R/a) 
g % hn s S he a \ ixe 

Here e 2 /hc ~ 1/137 is the fine structure constant. Physically, 1/g is the dimensionless 
zero-frequency impedance of the superconducting wire 

1 2uh 

In the absence of fluctuations, the phase varies linearly with the distance along the loop; 
therefore, it is convenient to introduce a new variable, x( x )> by the relation: 

x 

= (Pojr + Xix), (7) 

where (fo is determined such that the energy U has its minimum at \ = 0. This condition 
can be written in a simple form: 

L* L* 

siny^o + — y?o = -7-$. (8) 
For later use, we rewrite the expression @ for U in terms of \ as follows: 

U = Ej{- cos[x(27r) - X (0) + vo] - [x(2tt) - X (0)] sin cp 
L* in fdxV L 



Here a new ("angular") coordinate 6 = 2ttx/L has been introduced. 
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The DC Josephson effect at zero temperature is described fully by the ^-dependence of 
the ground state energy, E gT , of the system under consideration. Because only the term U 
of the energy depends explicitly on $, the persistent current </($) = (2e/h)dE sr /d$ can be 
expressed in terms of the average (U) over the ground-state wave function: 



III. QUALITATIVE ANALYSIS 

The behavior of J as a function $ is very different in the two limiting cases L <ti L* and 
L ^> L*, which we will discuss qualitatively below. 



A. The case L < L* 



For a relatively short loop, L < i*, the kinetic energy of the supercurrent dominates 
over the Josephson energy, and the phase difference (fo across the junction is completely 
determined by the flux threading the loop, (fo ~ $, as can be seen from Eq. (||). Classically, 
i.e., in the absence of phase fluctuations (x = 0), the dependence of the persistent current 
on $ is given by = J,?sin$, with J° = 2eEj/h. We will see below that this classical 
result holds in the limit 1/g — > 0. 

In the quantum case \ 7^ 0, we can neglect the effect of the Josephson junction on phase 
fluctuations when calculating the average in fllPP for L <C L* . In other words, for a short 
loop only the term with the integral in (|9|), which corresponds to the kinetic energy of the 
supercurrent along the loop, is important. Substituting Eq. (|9|) into Eq. fllPP, and neglecting 
terms G(L/L*), we find 



J($) = J c °(cos(x(27r) - X (0))> sin$. 



(11) 



The evaluation of the average (cos(x(27r) — x(0))) does not differ in fact from the well-known 
calculation of the Debye- Waller factor |2"3|] . We quantize the fluctuating field m the 
standard way as follows: 



1^1 h& 



(12) 



where operators a n satisfy canonical commutation relations. Using (|T2|), it is straightforward 
to evaluate (cos(x(27r) — x(0))) where the average is taken with respect to the quadratic 
Hamiltonian 



n=l 



hvpinn 
L 



a\a n + - 



At low temperatures and L 3> hv p \/A, we finally obtain </($) = J c sin $, with a renormalized 
critical current 



hv 



pi 



1/9 



LA J 



sinh (nLkBT/hvp\) 



1/9 



(13) 
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The energy A in fllTf ) appears as a high energy cutoff for the plasmon waves; the result (p!3|) 
holds for temperatures ksT <C A. The classical result J c — > J° is recovered in the limit 
l/g^O. 

At this point we would like to note that the Lagrangian defined by Eqs. (|]) and @ has 
been derived in the limit of small phase fluctuations, ((V%) 2 ) <C l/£ 2 (0). This poses an 
upper bound on the allowed values of 1/g, 

llg<< f (kFaf ^m. (14) 

The right hand side of this inequality is proportional to the number (/c^a) 2 of quantum 
channels in the wire, and typically is large. We will be interested in superconducting wires 
characterized by g ~ 1. For such wires, (|14|) is satisfied, and the local quantum fluctuations 



of the current due to the propagating plasma mode [12] are smaller than the critical current. 



We therefore neglect phase slip events p4| , p5f in the wire 



Let us estimate the exponent 1/g for an Al wire. If the wire is very dirty, with a mean 
free path I ~ 1 nm, we estimate the zero temperature coherence length to be £(0) ~ 40 nm 
and the London penetration depth to be Ai(0) ~ 500 nm. Present-day technology enables 
one to fabricate wires with a cross-sectional area S = a 2 ~ (50 nm) 2 . Such wires would be 
characterized by an exponent 1/g ~ 1.2. For cleaner wires with / ~ a, the exponent can be 
expressed as 1/g = (16/a(nm)) 3 / 2 . We conclude that typical values of the exponent should 
be in the range 1/g ^ 1. 

Quantum fluctuations suppress the maximum Josephson current below its mean-field 
value J°. This suppression depends on the loop length L; according to (|13|) J c — > when 
L —>■ oo. As we will discuss below, this is an artefact of the lowest order of perturbation 
theory, where the effect of the junction on the fluctuations in the attached wire is disregarded 
completely. 



B. The case L > L* 

As we have seen above, if L L*, the solution <po(§) of equation (|]) that provides the 
absolute minimum of energy varies continuously with $. For a large loop with L ^> L*, the 
solution (po(§) has discontinuities: 

if < $ < 7T, 

if 7t < $ < 2vr. (15) 

Because the Josephson energy dominates over the kinetic energy of the supercurrents, the 
phase <£>o remains "pinned" to the minima of the cosine potential. Correspondingly, in the 
absence of fluctuations the equilibrium persistent current J($) has cusps, 

J($) ~ J c °^r$ if < $ < 7T, 
J($) ~ J c °^r($ - 2tt) if 7T < $ < 2tt. (16) 



L* 

L* 

if ~ 2tt + - 27r) 

Lj 
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We expect quantum fluctuations (i) to renormalize the bare Josephson energy and thus to 
suppress the slope of the saw-tooth dependence; (ii) to smear the cusps at $ = (2n + l)n, 
as quantum tunneling will remove the degeneracy between pairs of states having the same 
values of energy but different values of ipo. This is shown schematically in Fig. [I]. 

However, due to the fact that the Josephson energy is not a weak perturbation if L > L*, 
we have to take its effect on the quantum fluctuations into account. This is very similar 
to the problem of pinning of a ID crystal [jnj. The decrease in (cos[x(27r) — x(0)]) with a 
growing length of the loop L should saturate when L exceeds the characteristic length L*. 
The saturation occurs because the Josephson coupling pins the low-frequency modes, thus 
preventing the logarithmic divergence of the phase fluctuations at the junction. 



IV. RENORMALIZATION OF THE JOSEPHSON ENERGY 

In this Section we will analyze the renormalization of the Josephson energy by the charge 
fluctuations in the framework of the renormalization group (RG) approach. This approach 
will enable us to treat both cases L < L* and L > L* in a unifying manner. We will 
restrict our analysis first to zero applied flux, $ = 0. It is convenient to perform a Wick 
rotation to imaginary time r and to consider the euclidean action S for the fluctuating 
field x{0, t) = x(0, t) — x(27r — 9, r), which can be easily obtained from the Lagrangian C 
Next one integrates out the fluctuations in r) away from the junction, and obtains the 
effective action for the field x(9 = 0), 



hq f duj 

4-7T 



j^|xM| 2 M- JdTE° jC osx(9 = 0,T). (17) 



This action can be studied by a standard perturbative RG method fDjfl . We introduce a 



running cut-off energy /z, and find a flow equation [26] for the dimensionless Josephson 
coupling energy Ej = Ejj\i: 

^ = (l-l/g)Ej + 0(E 3 j), dl = -dn/n. (18) 

This equation describes how the Josephson coupling Ej is renormalized when high-energy 
degrees of freedom are integrated out [p7P - From Eq. (|1|) it follows that, upon decreasing 



/i, the energy Ej flows to zero if g < 1, whereas Ej increases if g > 1 [pS j. Note that in the 
latter case perturbation theory breaks down as soon as Ej ~ 1. 

In order to investigate these cases in more detail, we integrate Eq. ([18|) from the high 
energy cutoff /i/j = A at which Ej = Ej down to a value fii = %v v \/Iq, characterized by some 
length l Q . As a result, we find 



1/9 
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The case g < 1. In this case, the result ([19]) remains valid for \xi — > 0, i.e., for the largest 
possible values of l . Putting l ~ L we thus recover our earlier result (0) at T = 0. At any 
finite length L > hv p \/A, the renormalized Josephson energy Ej is smaller than the plasmon 
level spacing of the loop, %v p \/ L. Therefore the perturbative analysis of Section [III A| applies, 
and we find for the flux-dependent Josephson current 

We see that quantum fluctuations will completely suppress the Josephson current as L — > oo 
i£g<l. 

The case g > 1. In this case, the situation is quite different. The RG procedure should 
be stopped when Ej ~ 1, i.e., when the cutoff energy /i reaches a value /xj which satisfies 
the condition m = Ej(pi). As a result is renormalized down to a value Ej S , and should 
be determined from the condition of self-consistency, 

E ?=^m ' (21) 



which yields 



Ef = E«(^) . (22) 



The value Ef is reached for l ~ (L*/p)(A/E3) 1/(ff_1) . We thus conclude that if 5- > 1, 
the result ([13]) at T = holds as long as L ^ L*; for larger values of L, the decrease of 
the Josephson coupling slows down, and eventually Ej saturates at the value E*f given by 
Eq. (|22"p. Further suppression of the Josephson energy is prevented by the fact that the 
modes x( a; ) & t frequencies u < Ej S /h are pinned by the Josephson coupling, and hence 
cannot participate in the renormalization. 

When using A as an upper energy cut-off in the derivation of Eqs. (|T9|) - (0), we assumed 
the condition (ffl) to be satisfied. In fact, the above treatment remains valid, even if (ffl) is 
violated, but the weaker condition L* 3> C, holds; in this case, A should be replaced with 
E c in (ID - (|22|). 



The result ( fig) for the flux- dependent Josephson current in a loop with L > L* remains 
valid for values of flux away from the cusp at $ = 7r; we just should replace Ej by Ej. The 
behavior of J(3>) for $ ~ 7r is strongly affected by quantum tunneling, which we will study 
in the next Section. 



V. QUANTUM TUNNELING OF PHASE 

As we have seen in Section [llllj| , in the classical limit the flux dependence of the persistent 
current has cusps at $ = (2n + l)n for a large loop, L ^> L*. At these values of $, a 
degeneracy occurs: Two states having a phase difference at the junction given by ip = 2mr 
and tpo = 2{n + 1)tc respectively, have the same energy. This degeneracy may be lifted by 
the quantum fluctuations of phase at l/g ^ 0. Tunneling between the two macroscopic 
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states characterized by different values of ipo induces a shift 8E of the ground state energy 
of the system. As a result, the cusps in the function J($) will be smeared. We will show 
below that the tunnel splitting 5E <C hv p \/L, i.e., it is smaller than the gap between the 
degenerate ground state and the first excited state of the loop with the junction. Thus, at 
zero temperature, we are dealing with an effective two-state system, and the flux dependence 
of J(5Q = $ — 7r) near <3> = ir will be given by |29[ 



J(M>) » 

tl Lj 

xil- . = }■ (23) 

In particular, we see that the smearing is characterized by a width 5& s ~ (L/L*)(5E/Ej), 
see Fig. [p. 

The tunnel splitting SE is proportional to the amplitude t for tunneling through a barrier 
of height Ej. This tunneling involves a varying phase field ip(x) along the loop, and therefore 
occurs in a mult i- dimensional potential landscape. The dominant contribution to the tun- 
neling action grows logarithmically with the system size L. For large L the amplitude t can 
thus be obtained within the WKB-approximation. However the pre-exponential factor has 
to be retained, as the leading term of the WKB-approximation yields a power-law, rather 
than an exponential, decay of t with L. The calculation of 5E is therefore conveniently 
performed with the use of instanton techniques [T^| which generalize the WKB-method to 



higher dimensions and enable one to evaluate the pre-exponential factor directly. 

We start our anaysis by introducing a new phase variable <p(x) = <f(x) — (p(L — x) + 
(2x/L — 1)$ (note that <p(L/2) = 0). In imaginary time, the Lagrangian for this phase field 
reads 



L/2 

C=^ I dx 
Air 



fp! 1 dr j 1 dx 



+£/cos0(O,r). (24) 



Here, we put $ = it, i.e., we consider the degeneracy point. The Josephson coupling energy 
Ej appearing in (E5) is assumed to be renormalized by the high energy degrees of freedom of 
the loop. It is a complicated problem to actually perform such a renormalization procedure, 
due to the strong anharmonicity of the potential U (Eq. (0)) close to the degeneracy point. 
For our subsequent treatment of the quantum tunneling process the following qualitative 
description of the renormalization scheme will be sufficient. We imagine to integrate out 
high-energy degrees of freedom in the spirit of Section 0, starting from A down to a cut- 
off [Ai, which satisfies the inequality Ej S <C [i\ A. These degrees of freedom are so fast 
that they "follow" the tunneling process and merely adiabatically renormalize the Josephson 
energy Ej down to the value Ej, where Ej S < Ej < Ej. As will be discussed below, the 
tunneling process itself consists of a slow part and a fast part. The former involves the 
remaining low energy degrees of freedom with energies up to Ej, whereas the latter involves 
those with energies between Ej and jii. 

We will be interested in tunneling between the initial phase configuration <pi(x, — T/2) = 
— (2tt/L)(x — L/2) at time r = —T/2 and the final configuration (f)f(x, T/2) = (27r/L)(x — 
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L/2) at time r = T/2. Both are classical configurations which minimize the energy U, see 
Eq. (H). The tunneling amplitude is characterized by the matrix element 

{^^e-^I^^N Jv<Pe- s ' h . (25) 

Here, AT is a normalization constant, H the Hamiltonian, and S = JL-f/ 2 drC(r) the action. 
The matrix element fl25| ) can be used to determine the ground state energy E'gr, because it 
decays as e~ E & r l n for T — > oo. 



The matrix element (|25|) will be evaluated in the so-called dilute instanton gas approxi- 
mation |I6| . We first will construct a single instanton (SI), i.e., a classical trajectory in the 
inverted potential — U between the configurations fa and 0/ that passes once, at a time r c , 
through the minimum of —U. According to Refs. [HJI5], the tunneling process consists of 



a fast and a slow part. The actual tunneling at the junction (i.e., the passage through the 
minimum at r = r c ) happens within a short time ro, and involves a part of the phase field 
with a length Lq = v p \Tq. The length Lq is determined by minimizing the total action of 
the SI; as we will see below Lq ~ L*, in agreement with [p4| , [15| . The rest of the phase field 
makes the transition fa — > 0y slowly. Therefore the SI consists of three steps, see Fig. |2|: 

(i) —T/2 < r < t c — r /2: slow adjustment of the phase field away from the junction 
from the initial configuration fa to the intermediate configuration <p(x > Lq) = 0; 

(ii) r c — To/2 < r < r c + r /2: fast tunneling at the junction involving a part of the phase 
field with length L ; 

(iii) t c + t /2 < r < T/2: slow adjustment of the phase field away from the junction 
from the intermediate configuration to the final configuration <pf. 

The matrix element (|25]) for a SI can be written as a product of two amplitudes, one 
corresponding to the slow and one corresponding to the fast contribution. 

We describe the slow adjustment by decomposing the phase field 0(x, r) into modes on 
the loop, 

4>(x,t) = 4> n ( T ) sin (2nnx/L). (26) 

n=l 

The upper cut-off n max ~ L/L indicates that slow adjustment involves the phase field away 
from the junction (x > Lq) only. For the initial (final) state of each mode we have 0„(r = 
— (+)T/2) = — (+)2/n; the dynamics of the modes is determined by the Lagrangian ( |24|) 
without the cos-term. This is a quadratic problem and the contribution to the matrix 
element (E5l) can be calculated exactly. We find 



™ max im,u„ 



n=l 



(h\e- HT/h \<P,)l, = n A/^^e-^, (27) 



where m s = ITigL j [txv v \) is the "mass" of each mode and uo n = 27rnv p \/L its frequency. The 
SI action for slow adjustment of the phase is easily calculated to be 

SfL~hglog(L/L ). (28) 

The fast tunneling of the phase field close to the junction (x < Lq) can be described by 
the following Ansatz fl4| : 
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(p(x, r) = o (t)[1 - x/Lq] ; - ir < (fi < n . 



(29) 



Substituting this into fl24|) , we find the Lagrangian for 0o( r ) : 

/ 7 , \ 2 



Air 



dr 



+ COS I 



(30) 



This Lagrangian describes the "rigid" tunneling of a part of the phase field with length L 
in terms of the motion of a particle with "mass" m,f — HgL / (67rv v i) and "coordinate" (po 
in an inverted double-well potential V(<po) — ~ % w pi0o/(47r-Lo) — Ejcos((p ). The action 
corresponding to a single passage through the minimum of V at <p = can be estimated to 
be 13 



Silt - EjL /v p i + ahg, 



(31) 



where a is a constant of order unity. Minimizing the total SI action Sf 1 = Sf{ ow + Sf^ t 
with respect to L we find L = hgv p \/Ej ~ L*. Following the standard treatment for 
tunneling in a ID double well potential outlined in Ref. one can easily estimate the fast 
SI contribution to (|25|), 



-HT/h\ 



^)fast 



2L* J 



exp 



a SI 

O- 



fast 



(32) 



The total matrix element ( |25| ) can now be calculated in the dilute instanton approxi- 
mation | 16| . One sums over all configurations of single instantons and anti-instantons that 
involve transitions (pi — > (pf and (pf — > (pi respectively. This is done under the assumption 
that Si's do not overlap, which is justified as long as Sf 1 3> h, i.e., for (L/L*) 9 3> 1. As a 
result we find 



(<Pf\e~ HT/n \&) 



2L* 



n 

n=l 



7rh 



x exp 

From the behavior of 



VpiT 



L* 



exp 



L* e 



(33) 



for T — > oo we infer that there are two low-lying energy 



eigenstates with energies E^. = E^ r ± 5E, where E® r is an irrelevant reference energy and 



5E 



gr 

hv 



pi -sf 



hv 



pi 



T 



(34) 



We see that indeed SE <C hv p \/L for large L, consistent with the assumption leading to 
Eq. (H). 
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VI. DISCUSSION 



The interplay between disorder and Coulomb correlations strongly influences the proper- 
ties of low- dimensional superconductors. This is well-known for thin films |3(| : A transition 
from superconducting (S) to insulating (I) behavior occurs upon decreasing the thickness of 
the film. Only recently, developments in fabrication techniques made experimental studies 
on in situ grown quasi- ID wires [31,32] possible. These indicate that a similar transition 
might occur in a superconducting wire upon decreasing its cross-sectional area S. Cor- 
respondingly, one may expect that if the parameter g is smaller than a certain threshold 
value g c , the wire should behave as an insulator on length scales even shorter than the loop 
circumference L, and the theory presented above ceases to be valid. 

Recent attempts to extend the description of T c — suppression in homogeneous thin 



films [30] to include quasi-lD homogeneous wires have met with considerable difficulties [33 



On the other hand, in a ID boson system disorder is known to induce a localized-delocalized 
transition, which occurs for strongly attractive interactions between the bosons [|34]]. More 
specifically, in terms of our model, the results of [34| would correspond to a transition to 
insulating behavior at a value g c = 3/2. This threshold in interaction strength is reduced 
in the case of two coupled chains H35| , and one may conjecture a reduction to g c — 1 for a 
multi-mode wire, in agreement with PS[ . 

An interesting model system which shows a S-I transition is a ID array of Josephson 
junctions pB| , |3T| . The behavior of such an array is determined by a competition of the 
Josephson coupling Ej between the islands (which favors a phase-coherent superconducting 
state) with the electrostatic energy Eq (which localizes Cooper pairs on the superconducting 
islands). If the array is superconducting, (Ej/SEo) 1 / 2 > 2ir, its low lying excitations are 



'phase waves' with a linear dispersion, oo(k) ~ yj8EjE k. We therefore speculate that a ID 
array containing a junction which is weakly coupled to its neighbors could be used to study 
the renormalization of the Josephson energy discussed in this paper. An advantage of this 
system is that the energies Ej and E , which depend on properties of the array, are well 
known and controllable in a typical experiment f38f ; in particular the Josephson coupling 
Ej can be chosen from a large range of values. This also would enable one to systematically 
probe the regime close to the S-I transition. 

In conclusion, we considered a thin superconducting loop which contains a Josephson 
junction. The "environment modes" of the loop, which consist of propagating plasmon 
modes with a soundlike dispersion, were found to renormalize the Josephson energy of the 
junction. The strength of the renormalization is determined by the dimensionless zero- 
frequency impedance of the loop. In order to observe this renormalization we propose to 
phase-bias the junction by threading the loop with a flux and measure the corresponding 
Josephson current. For a relatively short loop, the kinetic energy of the supercurrent domi- 
nates over the Josephson energy and the phase difference across the junction is completely 
determined by the flux threading the loop. The supercurrent depends on flux in a sinusoidal 
fashion. Quantum fluctuations suppress the amplitude of this dependence. In the opposite 
limit of a large loop, the phase difference remains more or less "pinned" to the minima of the 
Josephson energy. Correspondingly the persistent current shows a saw-tooth dependence on 
flux in the classical limit. Quantum fluctuations not only suppress the amplitude of the 
oscillations, but also affect their shape. While the impedance of the loop increases, the 
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cusps of the saw-tooth dependence are smeared; as the impedance tends to the quantum 
unit value the shape of the oscillations approaches a sinusoidal form. At any finite length 
of the loop, the renormalized current amplitude is finite and shows a power-law dependence 
on the junction conductance, with an exponent depending on the impedance of the loop. 
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FIGURES 



FIG. 1. Schematic dependence of the Josephson current J on $ for a large loop, L 3> L*: (a) 
sawtooth dependence, found in the absence of fluctuations; (b) quantum fluctuations suppress the 
slope of the sawtooth dependence and smear the cusps over a typical width e><3? s . 

FIG. 2. Four configurations of the phase field <p which occur during the tunneling process; these 
configurations are discussed in the text. 
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